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IMPEDANCE  OF  A  FINITE  INSULATED  CYLINDRICAL  ANTENNA  IN  A  COLD  PLASMA  WITH  A 
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by 

Janis  Galejs 
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Sylvania  Electronic  Systems 
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Waithem,  Massachusetts  0215  U 


ABSTRACT 

A  variational  formulation  is  developed  for  the  impedance  of  a  finite  cylind¬ 
rical  antenna  embedded  in  a  dielectric  cylinder,  which  is  surrounded  by  a  magneto¬ 
ionic  medium  (cold  electron  plasma)  with  the  static  magnetic  field  impressed  in 
a  direction  parallel  to  the  antenna  axis.  Closed  form  expressions  are  obtained 
in  the  limit  of  low  frequencies,  and  for  short  antennas  in  a  uniaxial  medium. 

The  impedance  of  a  short  antenna  is  nearly  the  same  as  for  an  assumed  triangular 
current  distribution,  except  that  further  resonances  are  observed  in  the  vicinity 
of  the  gyro  frequency,  where  the  antenna  becomes  electrically  long.  These 
resonances  may  be  shifted  to  frequencies  exceeding  the  gyro  frequency  in  the 
presence  of  an  insulating  layer  around  the  antenna.  For  very  thin  insulating 

layers  the  wave  number  of  the  variationally  approximated  current  distribution  is 

•  •  1  {  r-  '  1 

to  the  first  order  equal  kQ  ( is  the  leading  diagonal  element  of  the 

permittivity  matrix),  where  the  gyro  frequency  may  be  both  smaller  or  larger 
than  the  plasma  frequency.  However  this  approximation  does  not  apply  to  current 
distributions  along  the  insulated  antenna.  The  present  calculations  are  also 
compared  with  earlier  work  on  antenna  impedances. 
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Impedance  of  a  finite  insulated  cylindrical  antenna  in  a  cold  plasma  vlth  a 

longintudlnal  magnetic  field 

by 

Janis  Gale js 


1.  INTRODUCTION 

A  large  number  of  authors  have  been  concerned  with  radiation  from  various 

1-12 

sources  in  a  magnetoionic  medium.  The  impedance  of  a  short  current  element 

13-15 

has  been  considered  using  the  quasi-static  approximation,  and  integral 

expressions  for  the  radiation  resistance  of  elementary  antennas  have  been 

.  ,  .  16-2 6 
developed. 

The  impedance  of  a  linear  antenna  has  been  also  considered,  but  there  are 
numerous  discrepancies  between  the  results  of  various  investigators.  Ihe  imped¬ 
ance  of  the  cylindrical  antenna  with  an  arbitrary  directed  static  magnetic  field 

27 

has  been  computed  by  Ament  et  al.  for  frequencies  above  the  plasma  or  cyclotron 

frequencies.  This  Impedance  formulation  assumes  that  the  current  distribution 

28 

along  the  antenna  is  sinusoidal  vlth  an  assumed  wave  number.  Staras  has  computed 
the  impedance  of  a  short  antenna  (which  is  either  parallel  or  perpendicular  to  the 
static  magnetic  field)  based  on  the  concept  of  a  volumetric  current  distribution 
which  is  selected  such  that  the  evaluation  of  the  impedance  Integrals  is  reasonably 
practical.  Numerical  results  are  presented  only  for  frequencies  which  are  much 
lower  than  the  plasma  or  the  cyclotron  frequencies,  and  for  thin  antennas  parallel 
to  the  magnetic  field  the  antenna  reactance  is  inductive  and  it  is  1-1/2  to  8 
times  larger  than  the  resistance. 

Antenna  impedance  is  also  discussed  in  two  recent  reports.  Brandstatter 
29 

and  Penico  consider  an  arbitraxy  direction  of  the  static  magnetic  field,  and 
assume  that  the  antenna  current  can  be  represented  as  the  sum  of  two  sine  waves. 
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Th#  complex  vara  numbers  for  these  two  trial  functions  are  the  sane  as  for  the  plane 

waves  (ordinary  and  extra-ordinary)  propagating  along  the  antenna  axis.  Only 

numerical  results  for  frequencies  which  are  much  lower  than  the  plasma  or 

cyclotron  frequency  are  presented.  The  antenna  reactance  for  an  orientation 

parallel  to  the  static  magnetic  field  is  less  than  10  Ohms  for  antennas  of  the 

length  from  10  to  30  meters  at  the  frequency  of  10  Kc/s  and  the  resistance  part 

of  the  antenna  iuqaedance  is  3  to  U  orders  of  magnitude  lower.  This  seems  to  be 

28  30 

in  contradiction  with  the  results  of  Staras.  Ament  et  al.  also  consider  an 

arbitrary  direction  for  the  static  magnetic  field.  It  is  indicated  that  the 

antenna  impedance  is  capacitive  for  antennas  parallel  to  the  magnetic  field  in 

the  VI F  range.  For  the  antenna  parameters  investigated  by  Brandstatter  and 

Penico2%ie  antenna  impedance  exceeds  1000  Ohms  in  magnitude  (Table  2),  and  the 

impedance  has  a  significant  real  part.  Some  further  calculations  (Table  k)  show 

also  an  inductive  impedance.  Based  on  the  stationary  character  of  an  approxi- 

30 

mate  impedance  expression  Ament  et  al.  obtained  a  simple  estimate  for  a 
propagation  coefficient  of  the  current  distribution  along  the  antenna.  However 
in  the  low  frequency  limit  this  impedance  expression  differs  significantly  from 

15 

quasi -stationary  impedance  expressions  derived  by  others. 

The  discussion  above  indicates  that  further  work  is -needed  to  clarify  the 

impedance  characteristics  of  antennas  in  a  magnetoionic  medium.  The  impedance  of 

a  linear  flat  strip  antenna  in  the  presence  of  an  arbitrary  stratified  isotropic 

dielectric  has  been  computed  recently,  and  these  results  have  been  compared  with 

31 

experimental  data.  In  the  present  paper  a  basically  similar  variational  formu¬ 
lation  will  be  applied  to  a  cylindrical  antenna  geometry.  The  antenna  current  is 
assumed  to  flow  along  the  surface  of  a  cylindrical  shell  of  length  '2  I'  and 
radius  'a',  which  is  embedded  in  the  center  of  a  dielectric  region  of  radius  'b* 
within  a  magnetoionic  medium  with  tht?  static  magnetic  field  in  a  direction  parallel 
to  the  antenna  axis.  The  dielectric  cylinder  is  assumed  to  be  either  an  actual 


2 


dielectric  or  to  approximate  the  effects  of  an  ion  sheath  around  the  antenna. 

Such  a  region,  depleted  of  electrons  has  been  shown  to  affect  the  antenna  cnaracter- 

26 

istics  particularly  for  frequencies  below  the  plasma  frequency.  The  field 
exprefir  1'  in  the  anisotropic  medium  are  derived  using  one  dimensional  Fourier 
transforms.  The  considerations  are  restricted  to  a  slightly  lossy  medium  and 
the  radiation  condition  requires  only  that  the  fields  are  exponentially  attenu¬ 
ated  at  large  distances  from  the  sources.  The  fields  are  expressed  in  terms  of 
Bessel  functions  of  complex  argument  and  the  antenna  impedance  is  formulated 
after  satisfying  the  boundary  conditions  at  the  surface  of  the  current  filament 
and  at  the  interface  between  the  dielectric  and  the  plasma.  The  antenna  impedance 
is  computed  using  a  two-term  trial  function  for  the  antenna  current.  The  trial 
functions  are  assumed  to  be  two  sine  functions  which  will  be  specified  in  the 
later  discussion.  For  antennas  of  finite  radius  the  impedance  is  made  available 
in  numerical  form  only.  However  for  vanishingly  thin  antennas  closed  form 
expressions  are  derived,  in  the  low  frequency  limit,  for  short  antennas  in  a 
uniaxial  medium,  and  for  half  wave  antennas  in  free  space.  The  numerical  results 
are  compared  with  the  earlier  work  of  others  on  the  same  problem. 

2.  FIELD  EXPRESSIONS  IN  THE  COLD  ELECTRON  PLASMA  MEDIUM 

Assuming  a  (suppressed)  exp(-iut)  harmonic  time  dependence  of  the  fields  for 
a  source  free  medium  the  electric  and  magnetic  vectors  E  and  H  satisfy  the 
Maxwell's  equations 

V  x  E  =  iuno  H  (l) 

V  x  H  =  -  iue  e  •  E  (2) 

—  o  =  - 

where  u  and  €  are  the  permeability  and  permittivity  of  free  space  respectively, 
o  o 

For  a  z-directed  magnetic  field  the  dyadic  permittivity  is  of  the  form  tradition¬ 
ally  used  in  the  classical  magnetoionic  theory, 


.  4a*. 


€  = 


€1  iC2  ° 


”i€2  €1 


0 

e. 


where 


€,  = 


u>  2  (w+iv)/w 
1  +  -f 


=  - 


U)  2  -  (u*iv)2 

2  / 

U  U  /W 

P  c'. 

u,2  -  (l»H-iv)2 


w 


1  _  p 

u(w+iv) 


(3) 

(«0 

(5) 

(6) 


where  v  is  the  effective  collision  frequency  and  where  the  plasma  frequency 

and  the  cyclotron  frequency  u>c  of  electrons  are  defined  in  terms  of  the  electron 

density  N,  charge  e,  mass  m  and  the  applied  static  magnetic  induction  B  as 

o 

2  2 

w  =  e  N/(mc  )  and  u  =  |  e  |  B  /m.  Because  of  the  symmetry  of  the  excitation  the 

p  o  c  o 

field  components  will  be  independent  of  <t>  and  they  are  related  to  their  Fourier 
transforms  over  the  z  coordinate  by 

OO 

Fij(o'z)  =  -h  f  fij(o>w)  eivz  dv  (?) 

«-CO 

where  F  =  E  or  H,  i  =  p,*  or  z  and  where  subscripts  j  denote  particular  modes  of 
propagation.  A  substitution  of  (7)  into  (l)  and  (2)  gives  the  following  relations 


between  the  transformed  field  components 

1%  Hp  =  -  tv  (3) 

iUMo  “♦  =  -  4  h  +  iv  EP  (9) 

1(*o  \  ‘  S  ^  (o  \  >  (10) 

-iwH^  =  -  i<«0(ei  Ep  +  iea  Et)  (11) 

lw  «p  -  4  5z  *  *  iUe0  {'U2  Ep  +  61  (12) 

P  ^  (o  *♦)  -  ■  lu6o  '3  *. 


ll 


(13) 


vhere  the  arguments  (p,w)  and  the  subscripts  j  have  been  omitted.  Ihe  transverse 


field  components  of  (8),  (9).  (ll)  and  (l2)  can  be  expressed  in  terms  of  E  and  H 

Z  2 


as 


V  5 
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P 
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♦  iup 


where 


2,4 


k02  -  w2)  iW  55  \  +  %  c2  k02  1  «z  ] 

(14) 

2  v  |  -  v2)  ^  H,  ] 

(15) 

e2ko2  v  ^  -  lu>1o  <‘lko2  -  1,2  >  1  5z  ] 

(16) 

[-D  +  v2  (^/-v2)]  |j  \  *  lW0Z2ko2  ^  \  ]• 

(17) 

,  .  2  2\2 
-  (fjk,  -  v  ) 

(18) 

and  k  =»  €  .  Substitution  of  (l4)  to  ( 17)  into  (lO)  and  (13)  gives  two 


o  00 

coupled  wave  equations 


[-  D  +  v2  -  v2) 


$  E  +  a  2W  iun  i)  H  =  k  2  £-  D  E 
z  2  o  o  z  o  3  z 


6^k  2w£  E  +  lap  (€,k  2  -  w2)<$.H  =  iwp  D  H 

2  o  z  o'lo  '  z  o  2 

where  the  operator  Jb  is  defined  as 


(19) 

(20) 


*  0^  (p^) 

Substituting  H  of  (20)  in  (19)  H  is  related  to  E  by 
z  z  z 

iu%  'a  v  =  -1^  +  e3  KkoS  •  v2)  fz 


(21) 


(22) 


This  equation  for  H  can  be  differentiated  and  substituted  in  (19)  which  gives 

z 


the  following  fourth  order  differential  equation 


a 


-82e  +  b  Be  +  c  e  =0 


(23) 


where 


a  =  hV 


(24) 
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b  -  ko8  [koV  -  +  V35  -  "V  *3>] 


c  -  -  €3  kQ  D  (26) 

After  carrying  out  the  differentiations  indicated  by  the  operator  (2l),  the 

resulting  equation  can  be  exhibited  as  a  sum  of  a  fourth  order  differential 

equation  and  a  second  order  differential  equation  (equations  (4.26)  and  (2.162.1a) 

32\ 

respectively  of  Kamke  ) .Both  differential  equations  have  the  Hankel  functions  of 

/|i|\  2 

order  zero  Hqv  '(0p)  as  common  solutions  (m  =  1  or  2),  if  3  satisfies  the  quad¬ 


ratic  equation 
4  2 

a0  -  b0  +  c  =0 


(27) 


The  arguments  of  the  Hankel  functions  for  both  differential  equations  become 
identical  and  the  solutions  of  (23)  have  the  fonn 

E  =  B(v)  H  M  (0p)  (28) 

z  o 

This  result  can  be  verified  by  direct  substitution  of  (28)  in  (23)  when  noting 
that$H  ^(0p)  =  -  02  H  ^(0p).  The  solutions  of  (27)  are  given  by 


b  +  vb  -  4  ac 


For  w  =  0  the  above  solutions  0  correspond  to  purely  radial  wave  propagation  in 
a  direction  transverse  to  the  z-axis.  The  two  solutions  of  (29)  are  now  simplified 


1  (*!2  "  «22)  k02/*i 

g 

Following  a  convention  used  in  the  past  the  solution  with  a  plus  sign  in  (29)  which 
gives  0“  =  for  w  =  0  is  denoted  as  an  ordinary  wave  (subscript  +)  and  the 

other  solution  is  denoted  as  the  extraordinary  wave  (j=  -).  The  0=0  solutions 
of  (28)  represent  plane  waves  propagating  in  the  z -direction  with  wave  numbers 


w^  =  k  V  c,  +  €0 

+  O  x  —  d 
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In  the  presence  of  strong  magnetic  fields  «2  -*  0,  and  the  anisotropic  medium 
becomes  uniaxial  with  ^  f  «^.  Hie  two  solutions  of  (29)  become 


®+2  ■  ‘3(ko2  - 

b  2  =  a!  -  v2 
-  10 


(32) 

(33) 


Equations  (32)  and  (33)  apply  also  to  isotropic  media  when  «2<>  and  • 

The  ordinary  wave  (32)  corresponds  now  to  HI  field  components  which  satisfy  the 

wave  equation  ( 19)  •  Hie  extraordinary  wave  (33)  has  TE  field  components  and  it 

satisfies  the  wave  equation  (20).  Obviously  there  is  no  coupling  between  TE  and 

2 

TM  modes  when  e2  =  0.  For  a  lossy  plasma  3  will  be  generally  complex  and  the 

sign  of  3  can  be  determined  from  the  requirement  that  field  components  should 

decay  for  large  values  of  p.  Im(3)  should  be  positive  if  Hankel  functions  of 

the  first  kind  (m=l)  are  U6ed  for  the  solutions  (27)* 

In  an  anisotropic  medium  with  finite  magnetic  field  f*  0,  and  the  field 

components  E  .  and  H  for  a  given  mode  j  are  coupled.  The  coupling  it  character- 
ZJ  zj 


ized  by  the  admittance 


lJ 


H  « 

=  -5i 

E  , 

zj 


(3*) 


It  is  computed  after  substituting  (28)  in  (22)  as 

2  _  ,  2\  .  .  .2 


Y J  "  ‘  iwu 


i_  ci(fi.i  -  )  *  63vC 


€2V 


(35) 


Substituting  (29)  in  (35)  and  letting  e2  •*  0,  Y+  ~  *2  1*2  0  and  **z+  a  0 

Similarly  it  is  seen  that  |Y  I  ■*  «  and  E  =0 


due  to  an  excitation  by  E_ , 

z+  '  1  z- 

due  to  an  excitation  by  H  ,  as  can  be  also  anticipated  from  the  wave  equations 
(19)  and  (20). 


The  antenna  is  approximated  by  a  cylindrical  current  sheet  of  surface  den¬ 
sity  Jz  of  radius  "a"  and  total  length  "2i."  The  transform  of  the  current  density 

J  (p,z)  is  denoted  as  J  (p,w).  This  current  sheet  is  embedded  in  a  dielectric 
z  z 


7 
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cylinder  of  radius  ”b"  and  of  relative  dielectric  constant  .  The  dielectric 
cylinder  is  surrounded  by  an  electron  plasma  with  static  magnetic  field  impressed 
in  the  direction  of  the  cylinder  axis  (z -direction).  The  axial  field  components 


are  assumed  to  be  in  the  following  form 
Region  I  for  0  <  p  <  a 

\l  ’  B1  Jo<V>  (36) 

Region  II  for  a  <  p  <  b 

*  B2  [Ho<1)((5ie)  +  Vo(2)(V>]  <37) 

«Z2  •  Ai  Jo(eip)  <38) 

where  the  expression  of  remains  valid  also  in  Region  I  for  0  <  p  <  a. 

Region  III  for  b  <  p 

*  I  B3J  Ho(1)<V>  <»> 

J 

Bz3  -  Z  A3J  Ho(l)(ejp)  («» 

J 


b  2  2 

where  0i  -  w  and  0j  is  given  by  (29),  j  ■  +  or  -  ,  and  where  Jn(x) 

is  the  Bessel  function  of  the  first  kind  of  order  n. 

The  amplitudes  A^,  B^,  B^,  B^j  and  A^j  are  also  functions  of  w\  The  field 
components  E^  and  are  computed  after  substituting  (39)  and  (40)  in  (15)  and 
(17).  The  field  components  E^,  E^2  and  H^are  obtained  substituting  (38) 
in  (15),  (36)  and  (37)  in  (17)  and  by  noting  that  e2  =  0  and  e1  =  -  z±  in 

the  dielectric  region.  The  boundary  conditions,  H  -  H  =  J  and  continuous 

TU  fx  Z 

Ezl  =  Ez2at  p  =  a  can  com^ined  to  eliminate  B^.  After  employing  the  appro¬ 
priate  Wronskian  this  procedure  results  in 


J 

z 
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2  J^a) 


( “B2  +  h  *b> 


(41) 
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Eliminating  A.  from  the  equations  of  E  and  H  continuity  at  p  »  b 

1  f  3 


i  ■»[ 


YJ  -  *  blj  ~ 


*1 

r° 


Continuity  of  at  p  =  b  leads  to 


B2  [\<1)d»1«  +  W2)<^]  +I  B3jVi(1)(8jb) 


Continuity  of  at  p  =  b  gives 


'  B2  [Ho(l)(0ib>  +  Vo(2)<Bib>]  +I  -  0 


The  symbols  8,  b^  and  b^  are  defined  by 


4ic.k 

«  - 

iwp  n  0,  a 

o  l 

0*0*  r  ?  p 

b  =  — “rr  loop  (c,k  -  w^)  Y,  +  €_k 

lj  iwp  D  [_  o  1  o  '  j  2  I 


2v  ' 
o 

m 


b2j  *  ^  [  ko‘t(i2a  -  ‘l2)  +  6lko2''2  '  ^o^o8  W  YJ  ]  (47) 

with  D  defined  by  (l8).  The  four  equations  (4l)  to  (44)  contain  4  unknown 
amplitudes  B^,  BgRg,  B^+  and  B^_  ,  which  can  be  uniquely  related  to  the  trans¬ 
formed  current  density  J  .  The  antenna  impedance  will  be  specified  in  Section  3 
in  terms  of  the  axial  electric  field  E^  at  the  surface  of  the  cylindrical  current 
filament.  A  straightforward,  but  lengthy  algebraic  manipulation  shows  that 
E  at  p  3  a  is  related  to  J  by 

Zc  Z 


F(w)  =  — — • 


Q1  *o  "  Vl 
°ia4  ‘  a2a3 


where 


“l  *  W  -  W 


°2  =  al+a3-  '  al-a3+ 
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“3  ’  8  Ji(S1t')/J0(B1a)  (51) 

a4  *  8  Jo(P1l)/J0(eia)  (52) 

*o  ’  21  [  Jo<6la)  *o(f,ib)  *  V0ia)  Jo  <0lb)  ]  (53) 

\  *  21  [  J^*)  ^(8^)  -  lf0(81a)  Jx  (8jb)  1  (5b) 

H(l)(f»b)  H(l>(0b) 

alJ  ■  YJ  +  bU  -^bj-  <55) 

a2j  '  b2j  Hl(l)(0jb>  (56) 

a3j  ’  (£iko2/0i)  Ho<l)(f3.lb)  (57) 


and  6,  b,  b,..  are  defined  by  (45)  to  (47) .  Y  (x)  is  the  Bessel  function  of 
lj  2j  n 

second  kind  of  order  n.  F(v)  of  (48)  is  a  complicated  function  vhich  involves 


a  number  of  Bessel  functions.  For  uniaxial  anisotropy  (w  -♦  «)  vhen  e_  *  0,  (48) 
can  be  simplified  considerably.  It  was  indicated  earlier  that  Y+  *♦  0  and  Y  ~ 
if  €g  *♦. 0  .  Applying  these  conditions  and  noting  that  0+  is  given  by  (32),  (48) 


is  changed  to 


F(v)  = 


na  0i  Jo(0ia) 


U)  €  €. 

O  i 


«301  Hi(l)(0+b)  -  >!>,■,%)  \ 

-«l8+H0(l,(B+b)  JjSjb)  +  «3S1H1(1,(8+b)  JjSjb)  (58) 


F(v)  depends  only  on  the  wave  number  0+  of  the  ordinary  mode,  because  the  TE 

field  components  of  the  extraordinary  mode  are  not  excited  with  e2  =  0.  For  a 

finite  current  I  =  2na  J  finite  values  of  the  radius  *b*  and  a  vanishingly 
z  z 

small  radius  'a'  of  the  current  filament 


2na  J 


i0i2 

We  e. 
o  i 


(59) 


The  field  E  is  now  proportional  to  Y  (0  a)  which  becomes  logarithmically 

Z  0  1 

infinite  as  a  -*  0. 
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3.0  ANTF.Ni; A  IMPEDANCE 

The  driving  point  impedance  of  a  cylindrical  antenna  which  carries  a 
z-directed  current  may  be  computed  from  the  expression 


Z  =  - 


2na  /  E  J  dz 
•  z  z 

[I(z  =  0)]2 


(60) 


where  E  is  computed  at  the  surface  of  the  current  filament.  Z  of  (60)  is 

Z 

stationary  with  respect  to  small  changes  of  the  surface  current  density  J 

z 

or  of  the  current  l(z)  about  its  correct  value J  .  J  is  assumed  to  be  an  even 

z 

function  about  z  =  0.  Substituting  the  Fourier  integral  representation  (j)  for 
E  (z)  and  interchanging  the  order  of  integrations  it  follows  that 


Z  =  - 


[i(z=o) y 


a 


oo  oo 

/  dv  Ez (v)  /  J*A) 


cos  wz  dz 


(61) 


00  00 

J  dw  F(w)  J  Jz(z)  cos  wz  dz 


[iu=o)r 

•<«*  HI 

where  F(w)  is  defined  by  (48)  or  (58)*  The  antenna  current  density  is  assumed  to 
be  representable  by  the  trial  function 


J  (z)  =  A  sin 
z' 


kA(MzP 


+  B  sin 


kB(HzD 


(62) 


Substituting  (62)  in  (6l)  it  follows  that 


Z  = 


rAA  +  2  M  rAB  +  t.  7BB 
(AFa  +  BFb)2 


(63) 


where 


Fn  =  sin 

00 

7NM  =  '  T“  /  dw  F(w)  %(w)  ®m(w) 


itao 

l 


r  r  i 

%(w)  =  J  sin  ^(i-z)  cos  wz  dz  =  — - 2  (cos  w^-cos  k^g) 

O  kJJ  "  w 


(64) 

(65) 

(66) 


with  N,  M  =  A  or  B. 
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The  Impedance  (63)  is  a  function  of  the  amplitude  ratio  of  the  two  tried 


functions  (a/b).  Because  of  the  stationary  character  of  (60),  the  optimum 
ratio  (a/b)  is  determined  from  the  condition  dZ/d(A/B)  »  0.  This  leads  to 

A  7BBFA  "  7ABFB 


and 


B  7aafb  7abfa 


7aa7bb  "  7ab 


(67) 


(68) 


with 


’b  7AA  -  2FA  FB  rAB  +  FA  rBB 


(69) 


The  impedance  can  he  also  computed  for  a  single  term  trial  function  by  setting 
B  =  0  in  (63).  This  gives 


AA 
,  2 


(70) 


The  impedance  formulation  (68)  which  uses  the  antenna  trial  functions 

(62)  becomes  inaccurate  for  very  short  antennas  when  the  two  terms  of  the  current 

distribution  are  linearly  related  (i.e.  sin  kfiz  =  (k^/k^  sin  kAz  within  the 

computer  accuracy).  Inaccuracies  can  b6  also  anticipated  for  longer  antennas  if 

Fa  and  F_  approach  zero  simultaneously.  In  these  cases  sine  and  shifted  cosine 
31 

functions'^  have  been  used  (w  >  l.J  «  in  Figs.  3  to  5,  »  «  w  in  Fig.  8). 

P  p 


4.0  IMPEDANCE  OF  A  THIN  ANTENNA  WITHOUT  A  DIELECTRIC  LAYER 

For  an  antenna  without  an  insulating  shell,  b-*a,  ¥  of  (53)  is  zero  and 
^  of  (5*0  represents  a  Wronskian.  In  the  limit  of  a-*o  the  second  term  of  a 
predominates  and  it  follows  that  F(w)  of  (48)  can  be  expressed  as 


U 


H (l)(6  a) 

vM  -  a+  rnr 


-  d 


H, 


’(8  a) 


Hx(l,(e+a) 


(71) 
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where 


dj ' 


iwu  b 
o  lj 


bl+b2-  "  bl-b2+ 


(72) 


and  b^j  and  b^j  are  defined  by  (46)  and  (47). 


4.1  Quasistationary  approximations 

In  the  low  frequency  limit  k  -*  o,  and  it  will  be  assumed  to  be  negligible 

relative  to  the  integration  variable  w  in  all  of  the  previous  expressions.  The 

2  2  2  2  2 
wave  numbers  0^  simplify  to  0^  =  -  w  e^/e^  and  0  =  -  w  ,  which  makes  d+  =  o  and 

pmm  ■■  -  g 

d  =»  wuo w/(  v  )  .  The  antenna  impedance  is  computed  from  (70)  with  kA  of 

the  trial  functions  approaching  zero.  TMs  results  in 

m 

.2 


Z  = 
s 


2  j»2 

it  &Z 


f  dv  F(w/  ) 


(73) 


Substituting  the  above  values  of  d+  and  d  in  F(w)  of  (71)  and  applying  the  small 


argument  approximations  of  the  Hankel  functions 

00 


Z  =  - 

s 


2i 


2  2 
we  e,  it  Z 
o  1  o 


J  (l-cos  wi)2  log  w  +  log  (  |  ^ 


+  C 


(74) 


where  C  =  O.5772I  ...  is  Euler's  constant.  The  integral  can  be  evaluated  to  give 


z  = 

s 


we  e,  nZ 
o  1 


log 


( 1 ) 


1  +  \  log  Y- 
d  3  J 


(75) 


Z  of  (75)  is  the  same  as  the  quasi -stationary  approximation  (6 3)  of  Balmain.^ 
s 

It  will  have  a  resistive  component  even  for  a  lossless  medium  if  e^^  and  have 
different  signs. 


4.2  A  short  antenna  in  a  uniaxial  medium.  £k  Z  «  l) 

2 

In  a  uniaxial  medium  e1  =  1  and  €g  =  0,  and  the  substitution  of  0^  of 
(32)  and  (33)  in  (72)  shows  that  d+  =  0.  F(w)  of  (71)  is  simplified  to 


iwp  (k  -  w<L) 
_/  \  o'  o 

F(w)  =  — 


Ho(l)(0+a) 

HlU)(0+a) 


(76) 
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which  can  be  also  derived  from  (58)*  Substituting  (76)  in  the  integral  of  (70), 
assuming  kA  =  kQ  and  kQi  «  1,  and  using  the  small  argument  approximations  of  the 
Hankel  functions  the  real  part  of  the  impedance  Zg  is  computed  as 


r 


R  = 
s 


20  (k 


\  60rt/(koi) 


for  w  >  w 


for  u>  <  « 


(77) 


25 

This  expression  for  the  resistance  (77)  has  been  obtained  also  by  Seshadri.  J 


The  imaginary  part  of  Z  is  computed  as 

5 


log(f)  - 1  -  !iob  1*3!  ]  (78) 

which  is  valid  for  short  antennas  at  all  frequencies. 

The  quasistatic  approximation  (7C)  gives  the  same  R  and  X  values  as 

S  6 

(77)  and  (78)  for  to  <  This  similarity  between  the  quasistatic  approximation 

of  the  antenna  impedance  and  the  impedance  of  a  short  antenna  in  a  uniaxial 

33 

medium  has  been  discussed  recently  by  Mittra. 


W€  rti 


U.3  The  half-wave  antenna  in  free  space. 

2  2  2 

In  free  space  0+  =  k  -  w  ,  and  applying  the  small  argument  approxi- 

(1) 


mation  of  '(x),  (76)  simplifies  to 


F(v)  ■  ■  rS  f  k°2  ■ yS )  H°(1)  ( Jko2  • a ) 

o 


(79) 


For  a  half-wave  antenna  l  =  \/k  =  rt/(2  k  ),  k.  of  the  trial  functions  (62)  is 

O  n 

set  equal  to  kQ  and  the  ar^enna  impedance  is  computed  from  (70).  This  results  in 


Z  =  7 


AA 


.  2f  o 
00  dw  cos  (  2 

^  o 


.2  2 
k  -  w 
o 


•„“>  u 


,  2  2 
.,k  -  w  a 

M  o 


(80) 


In  the  limit  of  a  -♦  o  the  real  part  of  (80)  is  given  by 


Re  l  - 


(81) 
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which  is  equal  to  73.1  Ohms  following  equations  (36)  and  (37)  of  Galejs.  The 
imaginary  part  of  (80)  is  computed  as 


where  Kq(x)  is  the  modified  Bessel  function  of  second  kind  and  of  order  zero.  In 
the  limit  of  a  -*  o,  small  argument  approximation  of  the  Bessel  functions  can  be 
applied  to  both  integrals  and  (82)  simplifies  to 


This  is  equal  to  (-1+2.5)  Ohms  following  equation  (39)  of  Galejs.^1  Die  numerical 
values  of  (8l)  and  (83)  may  be  recognized  as  standard  results  of  antenna  theory. 

5.0  DISCUSSION  OF  NUMERICAL  RESULTS 
5.1  Uniaxial  medium 

The  magnetoionic  plasma  degenerates  into  a  uniaxial  medium  with  =  1 
and  =  0  in  the  presence  of  strong  magnetic  fields  (P  =  w  /u>  -*«>).  For  such 

c.  C  C  p 

a  medium  the  antenna  impedance  is  shown  in  Fig.  1,  where  the  left-hand  scale  of 

the  resistance  applies  for  w  <  u  and  the  right-hand  scale  for  w  >  w  .  The 

P  P 

calculations  are  made  for  antennas  of  length  X./20  and  \/l0  (indexes  1  and  2)  and 
for  very  small  (b  -a  <  10-1+  \)  and  for  larger  (b  -a  =  l/lOO)  radii  of  an  insulating 
layer  around  the  antenna  (indexes  p  and  i). 
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The  antenna  radius  a  and  length  2i  are  such  that  0=2  log  (2 i/a)  =8. 

The  medium  has  slight  losses  with  an  average  collision  frequency  v  =  0.001 

The  resistance  is  nearly  constant  for  w  >  increases  suddenly  near  w  =  and 

reaches  its  maximal  value  for  w  <  u^.  The  resistance  decreases  gradually  for  lower 

2 

values  of  u  or  higher  values  of  The  resistance  is  proportional  to  i  for 
u  >  Up  and  is  inversely  proportional  to  i  for  w  <  The  reactance  is  capacitive 
for  the  frequency  range  indicated  in  Fig.  2,  and  for  a  lossless  plasma  it  becomes 
logarithmically  infinite  at  u  =  wp.  The  reactance  is  decreasing  for  higher  values 
of  Up,  and  it  will  become  inductive  for  sufficiently  large  values  of  w^.  The 
same  general  impedance  behavior  is  also  observed  from  the  approximate  impedance 
expression  (77)  and  (78)  which  are  indicated  by  thin  lines  in  Fig.  1.  An  insulat¬ 
ing  layer  around  the  antenna  of  b  -a  =  \/l00  (index  i  in  Fig.  l)  has  little  effect 
for  frequencies  w  >  but  it  tends  to  decrease  the  resistance  and  to  make  the 
reactance  more  capacative  for  u  <  ■ 

5.2  Antenna  resistance  in  a  magnetoionic  medium. 


The  antenna  impedance  can  be  calculated  from  (70)  for  an  assumed  sinu¬ 
soidal  current  distribution  which  becomes  triangular  for  short  antennas.  The 
real  part  of  the  computed  impedance  is  shown  in  Fig.  2  for  an  antenna  of  length 
£  =  c/ojp  (£  =  l/(2n)  at  w  =  w^).  The  antenna  has  a  finite  radius  (ft  =  12.5)  and 
the  surrounding  medium  has  also  slight  losses  (  v/w^  =  10~ J) .  The  calcula¬ 
tions  are  made  for  a  plasma  of  ft  =  w  /u  =  0.5  and  1.5.  The  resistance  has  a 

c  c  p 

/  p  2~ 

large  but  finite  peak  at  the  upper  hybrid  resonance  frequency  wu  a*uc  +  wp 

which  is  indicated  in  Fig.  2  as  ft^  =  0^/10^.  The  radiation  resistance  becomes 

large  for  small  and  large  values  of  u/u^.  The  radiation  resistance  of  a  current 

filament  with  an  assumed  triangular  current  distribution  has  been  computed  recently 
22 

by  Seshadri  for  a  lossless  medium.  These  data  are  indicated  by  thin  lines  in 


Fig.  2 t  and  they  agree  with  the  present  computation  in  particular  for  the  lower 

frequencies.  However  for  a  lossless  medium  and  w  <  w  the  resistance  is 

c  p 
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equal  to  zero  in  the  frequency  range  w  <  w  <  u  ,  where  u  »  0.5  (-w  +  vw  +  ku>  ). 

C  li  J.  c  c  p 

The  resonance  peak  observed  for  w  =  1.5  w  at  w=<i)  is  higher  for  a  lossless  medium 

c  P  P 

and  the  two  sets  of  calculations  differ  significantly  for  higher  frequencies  where 

the  triangular  current  distribution  becomes  inaccurate.  However,  the  assumed  sinusoidal 
current  gives  a  high  resistance  for  l  «  x/2  (  or  near  u  =  3w^). 

5.3  Antenna  impedance 

The  antenna  impedance  is  computed  next  using  a  two  term  trial  function 

(62)  with  kA  =  kg/2  =  kQ  .  Such  trial  functions  can  be  shown  to  be  adequate  for 

impedance  calculations  in  free  space  if  the  antenna  length  is  not  near  a  multiple 

of  X/2,  and  they  will  also  indicate  the  frequency  ranges  where  the  electrical 

antenna  length  becomes  large  in  a  magnetoionic  medium  and  where  further  resonances 

may  be  expected.  The  calculations  of  Fig.  3  are  made  for  the  normalized  gyro 

frequency  fl  =w/w  =  0.5  and  an  additional  resonance  is  observed  for  frequencies 
c  c  p 

w  which  are  somewhat  less  than  u>c  ,  if  the  insulating  layer  is  of  negligible  thick- 

ness  (b  -a).  Increasing  losses(v/u>p  =10  )  will  decrease  this  resonance  peak  and 

the  losses  also  increase  the  resistance  in  the  frequency  range  w  <  w  <  w  ,  where 

the  resistance  is  zero  for  a  lossless  plasma.  An  insulating  layer  (or  ion  sheath) 

of  a  thickness  b  -a  =  X/lOO  shifts  the  resonance  to  frequencies  w  >  w  .  Also  it 

c 

makes  the  reactance  more  capacitive  for  w  <  and  slightly  more  inductive  for 

w  >  w  .  The  antenna  resistance  is  decreased  for  u  «  u  .  A  similar  impedance 
u  p 

behavior  is  also  seen  in  Fig.  4  for  fi  =  1.5.  The  additional  resonance  peak  is 

* 

observed  for  w  <  w  ,  if  the  insulating  layer  is  thin,  and  the  resonance  peak  is 

c 

shifted  closer  to  the  frequency  w  for  increased  thickness  of  the  insulating  layers. 

The  presence  of  the  resonance  peak  in  the  vicinity  of  00  indicates  that  the 

c 

antenna  becomes  electrically  long  as  w  approaches  w  .  Under  these  conditions  €, 

C 

approaches  infinity  and  the  resonance  can  be  explained  possibly  by  a  propagation 
coefficient  of  the  current  distribution  which  is  proportional  to  the  ,  follow - 

30 

ing  the  suggestion  of  Ament  et  al.  Further  calculations  are  therefore  made 


using  the  trial  functions  k  =  k  ,  k_  =  2  k  ,  where  k  =  vRe  e,  k  .  The 

A  s'  fB  s'  s  1  o 


1 
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Impedance  curves  shown  in  Fig.  5  exclude  the  frequency  range  w  <  w  <  w  ,  where 
Re  is  negative,  because  the  present  computer  program  is  intended  for  real 
coefficients  kA  and  kg  only  in  (64)  and  (66).  The  resonance  peaks  are  higher  than 
in  Fig.  3  and  4,  and  there  are  also  further  resonances  in  the  immediate  vicinity 
of  ii>c,  which  could  not  be  shown  in  the  frequency  scale  used  in  Fig.  5.  The 
impedance  computed  from  two -term  trial  functions  is  also  nearly  the  same  as  the 
computations  which  are  based  only  on  kA  =  kg  from  Eq.  (70).  The  ratio  of  the  trial 
function  amplitudes  A/B  is  computed  from  (67).  The  magnitude  of  a/b  is  of  the 
order  5  to  10,  and  this  shows  that  the  predominant  term  of  the  antenna  current 
representation  has  a  wave  number  k.  =  k  .  Further  calculations  have  been  made 

n  S 

using  k.  =  k  and  k  =  k  /2  .  Again  the  magnitude  of  A/B  is  larger  than  unity 
A  S  ij  6 

and  the  impedance  is  nearly  the  same  as  for  L  =  2  k  shown  in  Fig.  5*  This 

B  S 

suggests  that  the  wave  number  k  may  be  good  first  order  estimate  of  the  current 

s 

distribution  along  the  antenna,  if  the  thickness  of  the  insulating  cover  is 

negligible.  This  hypothesis  will  be  checked  by  plotting  several  variationally 

computed  current  distributions  and  by  comparing  them  with  a  single  sine  wave  having 

a  wave  number  k  . 

s 

The  calculations  of  Figs.  6a  and  b  are  made  for  frequencies  near  the 
resonance  peak  u>  <  u)c  .  For  small  values  of  (b  -a)  the  two  sets  of  calculations 
using  kA  =  kg  and  kfi  =  2  kg  (dashed  curve)  or  kg  =  0.5  kg  (dashed  and  dotted  curve) 
give  essentially  similar  current  distributions.  The  calculations  using  free  space 
wave  number  kA  =  kQ  and  kg  =  2  kQ  ,  indicated  by  solid  curves,  approximate  the 

* 

other  current  distributions  more  closely  in  Fig.  7b,  because  of  (4)  is  smaller 

in  magnitude  for  larger  values  of  u  ,  and  k  differs  less  from  k  for  u>  larrce. 

c  s  o  c 

The  thin  line  denotes  a  simple  sine  wave  with  a  wave  number  k  ,  and  it  lies  below 

s' 

the  other  curves.  This  indicates  that  the  wave  number  of  the  current  distribution 

is  somewhat  larger  than  k  in  Figs.  7&  and  7b. 

s 
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For  an  insulating  layer  of  thickness  (b-a)  =  l/lOO  the  current  distri¬ 
bution  is  quite  different.  It  is  nearly  triangular  for  wc/wp  =0.5  and  it  has  a 
relatively  low  peak  for  =  1.5.  Hie  insulating  layer  tends  to  make  the 

antenna  electrically  shorter  for  frequencies  near  and  the  simple  approximation 
of  the  vave  number  of  the  current  distribution  by  kg  that  may  be  Justified  for 
antennas  in  direct  contact  vith  a  surrounding  plasma,  definitely  will  be  in  error 
in  the  presence  of  an  insulating  layer  or  an  ion  sheath. 

Further  calculations  are  shown  in  Figs.  6c  and  d  for  antennas  of 
negligibly  thin  insulating  layers  and  for  frequencies  u)  in  the  immediate  vicinity 

of  u  .  Hie  thin  line  which  is  sine  wave  with  a  wave  number  k  ,  approximates  the 
c  s 

current  distributions  computed  with  =  kg  and  k^  =  2  kA  or 

agreement  is  particularly  close  for  uc/up  *  0»5>  but  even  for  wc/wp  =  1.5  the 

wave  number  k  gives  a  first  order  estimate  for  the  propagation  coefficient  of  the 
s 

current  distribution.  This  propagation  coefficient  k  should  not  be  identified 

s 

as  a  wave  number  of  surface  waves  which  do  not  exist  with  a  lossless  plasma  for 

34 

frequencies  w  above  the  lowest  of  or  u>c  following  Seshadri.  Obviously  the 

surface  waves  cannot  explain  the  resonances  shown  in  Figs.  3  and  5b,  which  occur 

in  the  frequency  range  <  u  <  uc  .  Also,  the  approximation  of  the  propagation 

coefficient  by  k  will  not  apply  to  antennas  insulated  from  the  surrounding  plasma, 
s 

It  has  been  suggested  that  the  wave  number  of  plane  waves  which  propagate 

29 

along  the  antenna  may  characterize  the  current  distribution  on  the -antenna. 

One  of  the  wave  numbers  given  by  (3l)  which  is  proportional  to  becomes 

very  large  as  w  approaches  wc  ,  it  is  nearly  real  for  u  <  u>c  and  it  becomes 
imaginary  for  u  >  w  .  However  for  w  «  u  its  magnitude  is  by  a  factor  of  >/sf 
larger  than  the  magnitude  of  k^ .  .Such  a  current  distribution  was  not  observed 
in  the  curves  of  Figs.  6c  and  d,  although  the  trial  functions  used  in  the  imped¬ 
ance  computation  (k  =k,k_=2k)  are  general  enough  to  account  for  it. 

A  Si)  S 


kg  =  kA/2.  The 
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5.U  Comparisons  with  other  impedance  calculations 

The  present  calculations  are  compared  in  Fig.  7  with  the  result  shown 

27  27 

in  *  *gs.  14  and  15  by  Ament  et  al.  for  w  >  w^.  Ament  et  al.  assume  a 

sinusoidal  current  distribution  with  a  wave  number  'Tel  k  .  The  two  sets  of 

3  o 

resistance  data  are  in  reasonable  agreement  but  the  differences  in  the  reactance 

curves  can  probably  be  explained  by  differences  of  the  trial  functions,  which 

should  have  the  most  pronounced  effects  near  the  resonance  length  of  the  antennas. 

The  present  calculations  are  further  compared  with  the  quasi-stationary 

approximation  ('/5)  using  the  plasma  parameters  shown  as  cases  2  and  3  in  table  4 
30 

of  Ament  et  al.  The  quasi-stationary  approximation,  indicated  by  thin  lines, 

follows  the  trend  of  the  more  accurate  calculations  and  it  shows  a  somewhat  lower 

antenna  resistance.  Hie  capacative  reactance  of  the  two  calculations  differ  by 

a  factor  of  nearly  two  in  the  case  B  .  The  data  shown  in  Fig.  8  differ  consider- 

30 

ably  from  the  calculations  of  Ament  et  al.  Their  resistance  values  are  in  most 

cases  lower  by  a  fc ’tor  of  approximately  100.  Furthermore  the  reactance  is 

inductive  and  is  three  times  larger  than  the  magnitude  shown  in  case  A,  and  it 

differs  by  several  orders  of  magnitude  for  most  of  the  points  of  case  B.  Hie 

30 

results  of  Ament  et  al  are  not  shown  in  Fig.  8. 

The  current  calculations  are  compared  further  with  the  quasi-stationary 

pQ 

approximations  and  with  the  computation  by  Staras  °  in  Fig.  9»  The  quasi- 
stationary  approximations  agree  closely  with  the  more  accurate  calculations  for 
long  thin  antennas.  However  for  short  antennas  (antenna  of  a  smaller  if  a.  ratio) 
the  quasi -stationary  approximation  indicates  a  resonance.  The  antenna  can  be 
shown  to  remain  electrically  short  in  the  vicinity  of  this  resonance  and  there¬ 
fore  the  accuracy  of  the  quasi-stationary  formula  should  be  questioned.  (This 
apparent  resonance  is  due  solely  to  the  decrease  of  the  if a  ratio  in  Eq.  (75) 
but  the  quasi-stationary  approximations  are  known  to  become  inaccurate  for 
antennas  of  decreasing  -C/a).  The  resistance  points  listed  by  Gtaras^  show  the 
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sane  trend  as  the  present  calculations,  hut  the  reactance  is  of  different  signs. 

28 

It  may  he  noted  that  Staras  uses  a  volumetric  current  distribution  vhich  does 
not  correspond  directly  to  a  fixed  antenna  surface  as  used  in  the  present 
investigation. 

29 

The  results  of  Brands  tatter  and  Fenico  were  briefly  discussed  in  Section  1. 
lhese  impedance  figures  are  of  such  magnitude  that  a  closer  c caparison  with  the 
results  of  the  present  calculations  does  not  appear  warranted. 

5.5  Limitations  of  the  analysis 

Only  antennas  of  1  <  0.51  are  investigated  in  the  present  paper.  For 

such  antenna  lengths  the  triangular  current  distributions  of  very  short  antennas 

are  inadequate,  but  it  is  not  necessary  to  apply  the  theory  of  long  antennas.  In 

this  intermediate  range  of  antenna  lengths  variational  Impedance  formulations  have 

been  successfully  applied  for  antennas  in  free  space  and  also  in  stratified 

32 

dielectric  layers. 

The  impedance  formulas  are  stationary  only  with  respect  to  small  changes 
AJZ  about  the  correct  current  density  Jcz,  and  large  changes  of  the  impedance  can 
be  anticipated  if  A  Jz  differs  significantly  from  Jcz .  The  converse  of  this 
property  can  be  readily  utilized.  If  the  Impedances  computed  using  different 
trial  functions  remain  nearly  constant,  it  may  be  concluded  that  the  assumed 
trial  functions  are  general  enough  and  that  the  resulting  current  distributions 
do  not  depart  significantly  from  the  true  current  density  J  ,  at  least  for 
purposes  of  calculating  the  impedance,  and  that  the  resulting  impedance  is  correct. 
The  constancy  of  the  impedance  can  be  verified  by  comparing  Figs.  3  and  4  with 
Fig.  3.  The  impedances  are  practically  the  same,  except  for  w  near  u  where  the 

-  k0  of 

Fig.  3  and  U  cannot  produce  a  sufficiently  oscillatory  distribution.  The  current 
distributions  of  the  variational  impedance  calculations  differ  by  approximately 
IQ  percent  in  Figs.  6a  and  b.  In  the  present  method  it  is  not  possible  to  select 
the  best  one  among  the  several  approximations,  although  an  auxiliary  condition 
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antennas  become  electrically  long  and  the  trial  functions  with  ■  k^/2 


of  S  "0  could  be  used  for  fields  on  a  perfectly  conducting  antenna. 

The  purpose  of  the  present  paper  was  to  show  that  the  impedance  of  the 

linear  antenna  can  be  computed  in  magnetoionic  media  with  current  distributions 

that  are  more  accurate  than  the  usually  assumed  triangular  or  half  ware  sinusoidal 

distributions,  and  to  emphasise  the  presence  of  multiple  antenna  resonance  as 

w  approaches  w  .  The  current  distributions  are  obviously  approximate,  but  the 
0 

computed  antenna  impedances  are  considered  reliable  vlthln  the  accuracy  of  the 
numerical  plots,  except  in  the  vicinity  of  the  frequencies  vhere  w  ■  w  .  Here 

V 

the  antenna  is  electrically  long  and  the  variational  Impedance  formulation  vill 
become  inaccurate  even  in  free  space  for  reasons  mentioned  in  the  concluding 
paragraph  of  Section  3. 

5.6  Summary  of  Conclusions 

This  paper  has  presented  Impedance  calculations  over  the  complete 

frequency  range.  For  a  uniaxial  medium  approximate  closed  form  expressions  compare 

with  the  more  accurate  numerical  calculations.  There  Is  good  agreement  with  the 

22 

antenna  resistance  as  computed  by  Seshadri  over  the  frequency  range  vhere  his 

triangular  current  distributions  is  Justified,  but  further  resonances  are  observed 

for  frequencies  in  the  vicinity  of  the  gyro  frequency,  fy  examining  a  number  of 

trial  functions  for  the  current  distributions  of  the  antenna  It  was  deduced  that 

represents  a  first  order  estimate  of  the  wave  number  for  the  current 

distributions  along  the  antenna  if  the  antenna  insulation  Is  of  a  negligible 

thickness.  This  estimate  applies  both  above  and  below  the  plasma  frequency  but 

it  is  not  valid  for  insulated  antennas.  For  frequencies  exceeding  the  plasma 

27 

frequency,  the  calculations  compare  with  the  work  of  Ament  et  al.  1  ,  and  for 
frequencies  much  lover  than  the  plasma  frequency  there  is  an  indication  of  the 

2ft 

validity  of  quasi-static  calculations.  The  resistance  values  listed  by  Staras 
can  be  compared  with  the  present  work.  However  he  indicates  an  inductive  react¬ 
ance  whereas  the  present  calculations  give  a  capacitive  reactance  for  the  same 
antenna  parameters.  There  is  a  definite  disagreement  between  the  present 
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29  V) 

calculations  and  the  recent  work  of  Brandstatter  andPenico  and  Anent  et  al. 
for  low  frequencies. 
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Figure  5.  Antenna  Impedance  Using  the  Trial  Functions  k4  =  k,/2  - 
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Figure  8.  Antenna  impedance  for  u  <<up.  Comparison  with  the 
quasistatic  approximation  (thin  lines). 
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Figure  9.  Antenna  Impedance  for  u  <<  Up.  Comparison  v/ith  the  quasistatic 
approximation  (thin  line)  and  with  computations  by  Staras. 
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